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We study a process of opinion formation in a population of agents whose interaction pattern is 
denned on the basis of randomly distributed groups of three agents, or triplets -in contrast to 
networks, which are defined on the basis of agent pairs. Results for the time needed to reach full 
consensus are compared between a triplet-based structure with a given number of triplets and a 
random network with the same number of triangles. The full-consensus time in the triplet structure 
is systematically lower than in the network. This discrepancy can be ascribed to differences in 
the shape of the probability distribution for the number of triplets and triangles per agent in each 
interaction pattern. 
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1. Introduction 

The representation of interaction patterns as graphs - 
or, as they are usually called, networks- is a widespread 
paradigm in science. Networks lie at the basis of the de- 
scription of a broad class of systems, ranging from intra- 
cellular molecular reactions and neural tissues, to artifi- 
cial objects such as the internet [Newman et ai, 2006; 
Boccaletti et ai, 2006]. These systems are conceived as 
populations of active agents (chemicals, neurons, com- 
puters), each of them occupying a network node. A link 
between two agents represents the possibility that they 
interact, mutually affecting their individual dynamics. 
Binary interactions, in fact, are an implicit assumption 
within the network picture. The inherent binary-like na- 
ture of all physical interaction laws may explain why, 
from the viewpoint of physicists, do networks provide 
such a satisfactory conceptual framework for applications 
of statistical physics to other branches of knowledge. 

It is not always possible, however, to reduce the dy- 
namical laws of a multi-agent system to binary inter- 
actions. Examples are found, for instance, in social dy- 
namics [Starkey et ai, 2000; Johnson, 2008]. Imagine a 
group of five people -A, E, I, O, and you- who meet to 
discuss and make a decision on a given issue. As soon 
as A exposes her views, your own opinion begins to re- 
act to them, perhaps changing a little, or being rein- 
forced. Most probably, however, before your reaction to 
A's views is definitively settled, you are already listening 
to E's. These have been influenced by A's speech and, 
in turn, are expected to modify both your preexisting 
opinion and the way it is affected by A's. As successive 
views are presented and overlap with each other, the in- 
tricacy of their mutual effect grows. The group's final 
decision will hardly be describable as the outcome of a 
sequence of binary events, each involving just two peo- 
ple. In other words, from the viewpoint of the group's 
dynamics, a representation in terms of a mere collection 
of binary links joining its members seems unsuitable. The 
intricate overlapping of the involved elementary processes 



suggests that the group should be considered as a single 
entity, subject to dynamical rules able to give an overall 
description of the relation between the group's states and 
its possible outcomes. 

In this paper, we consider a process of opinion formation 
in a population whose interaction pattern is defined in 
terms of randomly distributed groups of three agents, 
or triplets. Emphasis is put on comparison with the 
same process running on a random Erdos-Renyi network, 
where the role of triplets is replaced by triangles, each of 
them formed by three mutually linked nodes. The Erdos- 
Renyi network is built in such a way that the number of 
triangles equals, on the average, the number of triplets 
in the triplet-based structure. We find that the time 
needed to reach full consensus in this structure is sys- 
tematically shorter than in the network, in spite of the 
fact that the number of triplets per agent in the former is 
the same as the number of triangles per agents in the lat- 
ter. This leads us to compare in more detail the statistics 
of triplets or triangles in each interaction pattern. The 
analysis discloses the different effects of the two patterns 
on the dynamical process taking place in the population. 

2. Multiplet-based populations 

From the several possible ways in which a network can 
be characterized [Newman et ai, 2006], the one which 
is most suitable for the generalization we discuss in this 
paper is the traditional definition given by graph theory 
[Gross & Yellen, 1999]. A graph is defined as a set of 
vertices -or nodes, which we associate with the members 
or agents in a population- and a list of edges -or links, 
which represent the potential interaction between agents. 
Each edge is effectively defined by the pair of vertices 
that it connects. In other words, for a given population, 
a network is completely characterized by listing the pairs 
of nodes connected by links. Over the population, thus, 
a network is a collection node pairs. 
This characterization of the structure of a population 
can be immediately generalized. Instead of a collection 
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of node pairs (duplets) we specify a collection of node 
groups (multiplets), each group containing a certain num- 
ber m of nodes (m-plet). In this way, agents are aggre- 
gated into groups of different sizes. In networks, an agent 
may participate of several interaction pairs, depending on 
its connectivity. Similarly, in the mutiplet-based popula- 
tion, a given agent can belong to more that one multiplet. 
Multiplets are conceived to replace network links as the 
basic units associated with the elementary interactions 
between agents [Johnson, 2008]. 

Figure [T] gives two instances of small multiplet-based 
populations. In Fig. [T^, we have a 6-agent population 
whose interaction structure is defined by a duplet, (1,2), 
a triplet, (1, 3, 5), and a quintuplet, (2, 3, 4, 5, 6). Agents 
1, 2, 3, and 5 belong to more than one multiplet. Fig- 
ure [TJa shows a 6-agent population structured into three 
triplets, (1,2,3), (2,4,5), and (3,5,6). 
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FIG. 1: (a) A 6-agent population whose interaction pattern is 
given by a duplet, a triplet, and a quintuplet, (b) A 6-agent 
population, with three triplets, (c) A 6-node network, with 9 
links joining the nodes which share triplets in (b). 



Is there, however, a genuine difference between a 
multiplet-based structure and a network? Wouldn't each 
multiplet be equivalent to a fully connected group of ver- 
tices -a clique [Newman et al., 2006] or simplex [Johnson, 
2008]- in an ordinary network? Comparison of Figs. [TJa 
and c shows, by means of a simple example, that this is 
not the case. In Fig. QJ:, each triplet of Fig. [T]b has been 
replaced by the links which interconnect all its agents. 
In the resulting network, however, the triangle formed 
by agents 2, 3, and 5, turns out to be equivalent, from 
the viewpoint of its interaction pattern, to the triangle 
formed, for instance, by 1, 2, and 3. In the structure of 
Fig. Hh, on the other hand, (1, 2, 3) is indeed one of the 
triplets in the population, while (2, 3, 5) does not exist as 
such. 

A full appraisal of the implications of extending networks 
to multiplet-based structures must in any case be sought 
through its effects on the collective dynamics of the popu- 
lation. In the next section, we consider a process of opin- 
ion formation [Krapivsky & Redner, 2003; Boccaletti et 
al, 2006; Galam, 2008] and compare results for the time 
needed to reach full consensus in networks and in triplet- 
based structures. Differences can be ascribed to rather 
nontrivial statistical properties of both kinds of interac- 
tion pattern. 



3. Opinion dynamics in triplet structures 

The dynamical processes that may require a represen- 
tation of the population in terms of multiplets -for in- 
stance, the case of decision making discussed in the 
Introduction- need not be the same as those which work 
on networks -such as, for instance, peer-to-peer informa- 
tion transmission. For the sake of comparison between 
the effects that these two different interaction patterns 
have on dynamics, however, we analyze in this section 
the same process of opinion formation running both on 
a network and on a triplet-based structure. 
Consider a population of N agents where, at each time t, 
the opinion Si(t) of each agent i adopts one of two given 
values, say, Sj(i) = +1 or — 1. At each evolution step, 
three agents i, j, and k, are chosen, and they adopt the 
opinion of the majority among them, i. e. 

Si(t+1) = s 3 (t + l) = s fe (t + l) = 

= sign[ Sl (t) + Sj(t) + s k (t)}. (1) 

The choice of the three agents is performed according 
to the structure of the population, as explained in the 
following. 

Firstly, we consider an A-node Erdos-Renyi random net- 
work [Newman et al, 2006], where each of the N(N-l)/2 
possible links is effectively present with probability p. 
Thus, for large N, the expected number of links is 
L = pN 2 /2. In this network, we define a triangle as 
a set of three mutually connected nodes. The expected 
number of triangles is A = p 3 N 3 /6. On this network, the 
three agents chosen at each evolution step of the process 
of opinion formation must form a triangle. 
Secondly, we consider what could be called an Erdos- 
Renyi triplet-based structure, where each of the A 3 /6 
possible triplets in the A-agent population is effectively 
present with probability q. Consequently, the expected 
number of triplets is T — qN 3 /6. In the process of opin- 
ion formation running on this structure, the three agents 
chosen at each step must belong to the same triplet. Tak- 
ing q = p 3 , we have on the average T = A. With this 
choice, the dynamics on the network and on the triplet 
structure can be quantitatively compared. 
We study the process of opinion formation on these two 
interaction patterns by means of numerical simulations. 
We focus on the time needed to reach full consensus, 
when all the agents share the same opinion. Three stages 
contribute to the full-consensus time in a finite popula- 
tion. If in the initial condition the two opinions are more 
or less balanced over the population, the first stage corre- 
sponds to the escape from the vicinity of the state where 
opinions are equally frequent, to a region where the un- 
balance between opinions is of order N. In this stage, 
the dynamics is equivalent to a symmetric random walk. 
Therefore, its duration is of order A 2 . In the second 
stage -whose duration is of order A- the dynamics is 
essentially deterministic, and corresponds to the increas- 
ing prevalence of the majority opinion over the dissenters. 
Finally, when only a few dissenters remain, the dynam- 
ics is stochastically driven by the now-infrequent events 
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where dissenters are chosen and converted to the major- 
ity opinion. The duration of this final state is also of 
order N. In our numerical simulations, we skip the first 
stage by taking a number of dissenters equal to N/4, so 
that 75% of the population is already consensual. The 
initial dissenters are distributed at random. 
The numerical construction of the two interaction pat- 
terns is as follows. We first fix the number T of triplets, 
and choose at random the corresponding T groups of 
three agents. Two or more triplets containing the same 
three agents are avoided. With this procedure, the 
triplet-based structure is completely defined. Then, we 
construct a random network with L = iV(3T/4) 1 / 3 links. 
Multiple links between any two agents are forbidden. The 
expected number of triangles for that number of links is, 
precisely, T . For a given realization of the random net- 
work with L links, the number of triangles A -which is 
close, but not necessarily equal, to T - is determined by 
simple counting. The full-consensus time t c , measured 
in evolution steps, is determined as a function of A in 
the network, and of T in the triplet structure. Results 
for fixed networks and triplet structures, with given val- 
ues of A and T, are averaged over 10 4 realizations for 
different initial opinion distributions. 
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FIG. 2: Numerical results for the time needed to reach full 
consensus, t c , measured in evolution steps, for a process of 
opinion formation running on a population of N — 1000 
agents whose interaction pattern is a network (empty dots) or 
a triplet-based structure (full dots). The full-consensus time 
t c is plotted against the number of triangles A or triplets T 
in each pattern. Each dot stands for an average over 10 4 
realizations of the initial opinion distribution. 

Figure [2] shows t c as a function of A and T for a pop- 
ulation of N = 1000 agents. Empty and full dots re- 
spectively correspond to results for the network (t c vs. 
A) and the triplet structure [t c vs. T), for the number 
of triangles and triplets ranging from 2500 to 5000. We 
find that the number of evolution steps needed to reach 
consensus on the network is more disperse and system- 
atically larger than in the triplet-based structure. In the 



range displayed, the two values of t c differ by a factor 
of almost two. For larger values of A and T, as it may 
be expected, the full-consensus times for the two interac- 
tion patterns approach each other. In fact, in the limit 
where all the possible links and all the possible triplets 
are present, A = T = N 3 /6, the two times must be 
identical. Conversely, they are increasingly different as 
the interaction patterns become more sparse. Consistent 
results were obtained for systems of various sizes. 
What is the origin of the difference between the full- 
consensus times in the network and in the triplet struc- 
ture? In view of the fact that the dynamical rules of 
opinion formation are identical for the two populations, 
the answer is to be found in their interaction patterns. 
In the next section, we provide evidence that an excess of 
network nodes belonging to a small number of triangles 
with respect to the agents belonging to the same number 
of triplets explains the difference. This requires to ana- 
lyze in detail, for each kind of pattern, the distribution 
of triangles or triplets over the population. 

4. Statistics of triangles and triplets 

In our numerical simulations, we have compared the full- 
consensus times in a network and in a triplet structure 
constructed in such a way that the number A of triangles 
in the former is, on the average, the same as the number 
T of triplets in the latter. This implies, in particular, 
that the average number of triangles per agent, 3A/7V , 
equals the average number of triplets per agent, 3T/N . 
In the following we show that these two coincident aver- 
ages correspond however to different distributions. More 
specifically, the fraction of agents which belong to a cer- 
tain number of triangles in the network is in general dif- 
ferent to the fraction of agents which belong to a certain 
number of triplets in the triplet-based structure. 
The fraction of nodes which belong to exactly S triangles 
in an iV-node Erdos-Renyi network can be evaluated by 
first considering the probability P n that a generic node 
is connected to exactly n of the possible TV — 1 neighbors, 
which is given by the binomial distribution 



Pn = 



N - 1 



p n (l-p) 



N-l- 



(2) 



where p is the probability that a network link is effectively 
present (see Section 3). The n neighbors form n(n — 1)/2 
pairs, some of which can in turn be joined by links, thus 
forming triangles. The probability Ps(n) that the node 
with n neighbors belongs to 6 of those triangles reads 



Ps(n) 



n(n-l)/2\ s 



p A {l-p) 



n(n-l)/2— 5 



(3) 



The total probability Ps that, irrespectively of the num- 
ber of neighbors, a node belongs to exactly S triangles is 
given by the total contribution over the different values 
of n: 
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(4) 
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Some little algebra shows that, as advanced, the average 
number of triangles per agent is, in the limit of large N, 
(5) = p 3 N 2 /2 = 3A/N. The width of the distribution 
can be computed by approximating P n by a Gaussian, 
yielding for the variance 



-I = T P 3 (1 



p)[4N 2 p 2 -2N{p 2 +p -1)-1]. (5) 



Which term dominates this expression for large N de- 
pends on how p is assumed to depend on N in that limit, 
as discussed below. 

The probability distribution P T of the number of triplets 
per agent r in the iV- agent population with T triplets is 
found, more straightforwardly, to be the binomial 



Pr 
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N 



T-i 



(6) 



The mean value, as predicted, is (r) = 3T/JV. For large 
N, the variance is a 2 = 3T/JV. 
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FIG. 3: Distribution probability for the number of trian- 
gles per agent in a network (empty dots) and the number 
of triplets per agent in a triplet-based structure (full dots), 
measured numerically for a population of N = 1000 agents 
with 4500 triangles or triplets. Curves stand for the analyti- 
cal results of Eqs. @ and @ . 



Figure [3] shows the probability distributions Ps and P T 
for a population of N — 1000 agents with 4500 triangles 
or triplets. Curves stand for the analytical expressions of 
Eqs. ([3]) and ©, and dots correspond to the numerical 
evaluation of the probabilities, which turn out to be in 
full agreement with the analytical results. The average 
number of triangles or triplets per agent is the same for 
both distributions, {§) = (r) = 13.5. Regarding their 
shapes and widths, on the other hand, the two functions 
are substantially different. For the triplet structure, we 
recognize the symmetric binomial distribution of Eq. ^ . 
The distribution of triangles per agent, in contrast, is 
clearly asymmetric, with a rather long tail for large S. 



Since, moreover, the maximum of Ps is sensibly lower 
than that of P T , for small values of 5 and r the probability 
of triangles per agent is higher than that of triplets per 
agent. In other words, the fraction of agents which belong 
to a small number of triangles in the network is larger 
than the fraction of agents belonging to the same number 
of triplets in the triplet-based structure. 
The difference between P$ and P T for small values of 
their variables provides an explanation for our observa- 
tion in Section 3 that the full-consensus time in the net- 
work is systematically larger than in the triplet structure. 
Once the quasi-deterministic stage of the opinion forma- 
tion process has elapsed, only a few dissenters remain in 
the population. These last non-consensual agents are ex- 
pected to belong, on the average, to a small number of 
triangles or triplets. In fact, those agents belonging to 
a relatively large number of triangles or triples have had 
more opportunities of being selected during the preceding 
stage, thus reaching consensus earlier. For a given num- 
ber of agents in the left end of the distributions Ps and 
P T , however, the average number of triplets per agent is 
larger than the corresponding number of triangles. Con- 
sequently, the last dissenters have a larger probability 
per evolution step to be chosen, and their opinion thus 
changed, in the triplet-based structure. The final evolu- 
tion stage will therefore elapse faster than in the network. 
We end our discussion on the distributions Ps and P T by 
comparing their limits for large populations, N — > oo. 
As advanced above, to effectively perform the limit it is 
necessary to specify how the interaction pattern changes 
as N grows. We consider three cases: (i) constant to- 
tal number of triangles or triplets; (ii) constant average 
number of triangles or triplets per agent; and (iii) con- 
stant fraction of effectively present triangles or triplets. 
In case (i), A and T are independent of N, and the prob- 
ability p of the Erdds-Renyi network scales as p ~ N . 
In case (ii), both A and T must be proportional to N, 
while p ~ N~ 2 / 3 . In case (iii), A and T scale as iV 3 , and 
p remains constant. 

It is clear that, since in case (i) the number of triangles 
or triplets is kept constant as the population grows, the 
distributions Ps and P T in the limit of large N are triv- 
ial. They both collapse to delta functions, with vanishing 
mean value and variance. Case (ii) is more interesting. 
In this case, as N grows, the mean values (S) and (r), 
and the variances a 2 and a 2 approach constant values. 
By construction, the limiting values of (S) and (r) are 
identical. Also the two variances turn out to be iden- 
tical in the limit. Interestingly, however, their mutual 
approaching for large N is very slow, with 
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where A; is a constant. This suggests that the dynami- 
cal effects of the difference between Ps and P T -such as 
the different duration of the process of opinion forma- 
tion discussed in Section 3- may be perceivable even in 
very large populations. Finally, in case (iii) the limit- 
ing distributions are drastically different. In fact, for the 
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network the variance grows as c 2 . ~ A 5 , while for the 
triplet structure we find a\ ~ A 2 . 

5. Conclusion 

Multiplet-based structures are expected to replace net- 
works as a representation of the interaction pattern un- 
derlying a population when the relevant dynamical pro- 
cesses cannot be assimilated to sequential binary inter- 
actions. This paper, however, has been focused on a 
comparison of the dynamical effects of the two kinds of 
pattern on the same process of opinion formation, driven 
by a majority rule in groups of three agents. By means 
of numerical simulations, we have determined the time 
needed to reach full consensus, on one hand, in a pop- 
ulation structured on the basis of randomly distributed 
triplets and, on the other, on a population of the same 
size whose interactions are described by a random Erdos- 
Renyi network forming as many triangles as the number 
of triplets in the triplet-based structure. In spite of the 
strong similarity between the two versions of the opin- 
ion formation process, there is a noticeable difference in 
the full-consensus time measured on both interaction pat- 
terns. We have shown that the difference can be ascribed 
to higher-order statistical properties in the distribution 
of triangles or triplets per agent in each pattern. This 
conclusion points out the significant effect of the kind of 
underlying interaction pattern in the collective dynamics 
of the population. 

The present results should be considered as a preliminary 
step in the study of the dynamical effects of multiplct- 
based structures. Further work may address dynami- 
cal processes more specific to this kind of structure, for 
which networks do not provide a suitable representation 
of interactions. The possibility that in a given structure 
the size of multiplets in not homogeneous, but varies ac- 



cording to a prescribed distribution; the extension of the 
network-specific notion of degree distribution to multi- 
plets; the consideration of small-world or scale-free mul- 
tiplet structures are, among other issues, worth investi- 
gating. 
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